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Abstract 


The commonly used Hamiltonian of the chromomagnetic hyperhne splitting 
is inversely proportional to the product of the masses of two constituent quarks 
composing the meson. So it is expected to have (M^* — Md)/{Md* — ~ 

{Mb* — Mb)/{Mb* — MbJ — 1.6, when the constituent quark masses rUu^ = 0.33 
GeV and = 0.53 GeV are used. However, the experimental results show that 
the above ratios are very close to 1. We solve this puzzle by employing the Hamil¬ 
tonian recently proposed by Scora and Isgur. 
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Review of Particle Physics [Q] presents the masses of the D mesons (MeV): 

Md± = 1869.3 ± 0.5 , Mdo = 1864.5 ± 0.5 , M^± = 1968.5 ± 0.6 , 

Mz)*± = 2010.0 ± 0.5 , Mo*o = 2006.7 ±0.5, M^*± = 2112.4 ± 0.7, (1) 
Md*± - Md± = 140.64 ± 0.09 , Md^o - M^o = 142.12 ± 0.07, 

M^*± - Mjj± = 143.8 ± 0.4 , 

and the masses of the B mesons (MeV): 

Ms± = 5278.9 ± 1.8 , M^o = 5279.2 ± 1.8 , M^o = 5369.3 ± 2.0 , 

Ms* = 5324.8 ±1.8, M^jo = 5416.3 ± 3.3 , (2) 

Mb* -Mb = 45.7 ± 0.4 , Mb*o - M^o = 47.0 ± 2.6 , 

where Mb = {Mb± + Mbo)/2 . We note that Mj^*± — Mj^± is very close to M^*± — 
Mb,± and Mb,*o — Mbo in (|I|), and that Mb*o — Mbo to Mb* — Mb in (|^). 

The commonly used Hamiltonian for the chromomagnetic hyperhne splitting 
of vector and pseudoscalar mesons is given by 0 

where = ^^/dvr of the QCD running coupling constant g, and mq (sq) and 
rrig (Sg) are the masses (spins) of the heavy quark Q and the light antiquark q 
inside a heavy-light meson, respectively. The Hamiltonian in (^ corresponds to 
the potential due to one-gluon exchange for s-wave bound states. Treating the 
chromomagnetic hyperhne splitting as a perturbation, one obtains 

where 4 /q^( 0) is the two-body bound state wave function at origin. When it is 
assumed that Du,d and Dg (or Bu,d and Bg) have similar values of Ug and |4/Qq(0)|, 
Mm* — Mm is proprtional to l/mq [q = h, J, s) for D (or R) mesons. Then it 
is expected to have (Mb* — Mb)/{Mb* — MbJ — {Mb* — Mb)/{Mb* — MbJ — 
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1.6, when the constituent quark masses mu,d = 0.33 GeV and = 0.53 GeV are 
used. However, the experimental results in (2) and (|^) give the values 


M£)*± - M£)± 


Md*± 




= 0.978 ±0.003, 


Mq* — M, 


B 


= 0.972 ± 0.054, 


(5) 


Mb*o — Mgo 

which are very close to 1. Goity and Hon noticed this peculiar property of 
the heavy meson masses. Randall and Sather @] emphasized this discrepancy and 
called it a puzzle. They suggested that these data might be an interesting probe 
of heavy mesons. The purpose of this letter is to solve this puzzle by employing 
the Hamiltonian recently proposed by Scora and Isgur. 

Scora and Isgur 0 proposed the following Hhf(j) as a modification of Hhf{r) 
in 
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where the term in the parentheses would be the ordinary Fermi contact term in 
(H) if the anomalous coupling coefficient a were unity, and where Ei = {rrif ± p^)^. 
For s-wave bound states, the expectation value of the Hamiltonian Hhf{'r) in (^) 
is given by P] 
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where S' = 0 or 1 is the total spin of the meson. 

In order to obtain the heavy meson wave function $Qq(p) in (||), Scora and 
Isgur P] applied the variational method to the Hamiltonian with the nonrelativistic 
kinetic energy terms. They adopted the Gaussian wave function as a trial wave 
function for the heavy meson ground state. We follow the same procedure as 
theirs in obtaining the heavy meson ground state wave function, except for that 
we use the relativistic kinetic energy terms since the velocity of the light quark 
inside heavy meson is large. We apply the variational method to the relativistic 
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Hamiltonian 


H ^ + mg^ + l/(r), 


(9) 


where r and p are the relative coordinate and its conjugate momentum. We take 
the following potential energy for V{r) in (^, 

QJ _ 

V{r) = - 'i + Kr + Hhf{r), (10) 

r 

where Hhf{r) is the modihed chromomagnetic hyperhne Hamiltonian in (|^) pro¬ 
posed by Scora and Isgur. We take the Gaussian wave function as a variational 
wave function for the heavy meson ground state, 


vl;(r) = (^)3/2g-/3V/2^ 
'7T 


<h(p) = 


,-pV2/92 


The ground state wave function is then given by P, § 

dE{f3) 


{H) = {ij\H\ij) = E{(3), 


df3 


= 0 at /3 = /3, 


( 11 ) 


( 12 ) 


where (3 represents the inverse size of the meson = 3/(2 /3)), and E = E{]3) 

the meson mass. In explicit calculations, we used the following four different 
potential models which have the values of Oc, K and quark masses given in Table 
1: (A) Scora and Isgur P, (B) Eichten et al. 0], (C) Hagiwara et al. Q, and 
(D) Model D which is the same as the Hagiwara et alJs model (C) except for the 
values ac{D) = 0.48 and ac{B) = 0.32 which are given from the running coupling 
constants at the energy scales of Mr, and Mb, respectively. 

For each of the four models of (A), (B), (C), and (D), we obtained the function 
E{f3) in (P!2|). As a representative, we present E{/3) = M for the model (C) of 
Hagiwara et al. in Fig 1, where (a), (b), (c), and (d) correspond to Da, Dg, Bd, 
and Bg mesons, respectively. In Fig 1, the graphs with the mark □ are what we 
obtain when we take the anomalous coupling coefficient a in (^ as zero, that is, 
when we do not include the Hhf{r) term in (13). The values of the variational 
parameter fd which give the minimum of E{j3) are presented in the hld(O), T*s(0), 
Brf(O), and Bg(0) columns of Table 2. 
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The graphs with the marks O and + in Fig. 1 correspond to those for the vector 
and pseudoscalar mesons, respectively, which we obtain when we use desirable 
values of a, as we explain in the following. By including of the Hhf{r) term in (pH] ) 
with a fixed value of the anomalous coupling coefficient a in (^), we obtain E{(3) = 
M and its minimum from the condition in (|T^, for the vector and pseudoscalar 
mesons separately. Then, by performing the same procedure with different values 
of a, we obtain the vector and pseudoscalar meson masses Mm* and Mm which are 
given by the minimum of E{f3), and their difference AM = Mm* —Mm as functions 
of the coefficient a. We obtained these funtions for each of the four models of (A), 
(B), (C), and (D). In Fig. 2, as a representative, we present AM = Mm* — Mm 
as functions of a for Dd, Dg, and Bg mesons obtained for the model (C). 

The nice feature of Fig. 2 is that the graph of Dd is very close to that of 
Dg, also for Bd and Bg, and this interesting character is the clue for solving the 
puzzle in this letter. The dotted vertical lines indicate the values of a which give 
rise to the experimental values of AM. In Table 3, we present these values of a 
and AM = Mm* — Mm obtained by these a values. Table 3 shows that the ratio 
R = AM{Md) / AM [Mg) are obtained as the values which are very close to 1 for 
both D and B mesons, and the puzzle is solved. The graphs with the marks O 
and + in Fig. 1 are E{I3) = M obtained with the values of a in Table 3. Fig. 3 
shows that the minimum of E{f3) for the vector meson is bigger than that for the 
pseudoscalar meson as it should be, and the value of the corresponding f3 (denoted 
by dotted vertical lines) for the pseudoscalar meson is bigger than that for the 
vector meson. This means that the size of the pseudoscalar meson is smaller than 
that of the vector meson, since = 3/(2 /?). We present the sizes of the vector 

and pseudoscalar mesons and their ratios in Table 4. 

In conclusion, we have solved the puzzle of (Mp. — M^)/{Mr,* — M^J ~ 
{Mb* — Mb)/{Mb* — Mb^) ~ 1 by adopting the modified chromomagnetic hyper- 
fine Hamiltonian in (^ proposed by Scora and Isgur [Q instead of the commonly 
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used chromomagnetic hyperfine Hamiltonian in (j^). 
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Oic 

K 

rUu = THd 

TUs 

TUc 

rUb 

Scora 

0.67 

0.18 

330 

550 

1820 

5200 

Eichten 

0.52 

0.18 

330 

530 

1840 

5180 

Hagiwara 

0.47 

0.19 

330 

530 

1320 

4750 

Model D 

(T))0.48, (5)0.32 

0.19 

330 

530 

1320 

4750 


Table 1: The values of the parameters in the potential Eq. (IT) of the models used 


in this letter. 



DM 

D*d 

Dd 

DsiO) 

d: 

Ds 

Bd{0) 

B*d 

Bd 

BsiO) 

b: 

Bs 

Scora 

556 

526 

715 

592 

561 

767 

663 

646 

729 

716 

697 

794 

Eichten 

515 

488 

662 

547 

518 

709 

591 

577 

647 

636 

620 

701 

Hagiwara 

484 

460 

605 

512 

487 

642 

574 

561 

625 

616 

601 

673 

Model D 

486 

462 

608 

514 

490 

639 

523 

511 

567 

558 

545 

608 


Table 2: The values of the variational parameter f3 (MeV) of the Gaussian wave 
function which minimize {H). Here, the values in the column of i5d(0) are f3's 
obtained without taking Hhf{r) into account, and those in the columns of and 
Dd are ,d’s obtained with taking H^fir) into account for vector and pseudoscalar 
mesons, respectively, and the same for other mesons. 
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a{D) 

AM{Dd) 

AM{Ds) 

R(D) 

a{B) 

AMiBd) 

AM{Bs) 

R(B) 

Scora 

1.21 

141 

144 

0.98 

0.78 

46.0 

49.5 

0.93 

Eichten 

1.85 

143 

144 

0.99 

1.27 

45.7 

48.2 

0.95 

Hagiwara 

1.86 

147 

143 

1.03 

1.37 

45.3 

47.2 

0.96 

Model D 

1.80 

146 

143 

1.02 

2.49 

45.6 

46.5 

0.98 

Expts. 

— 

141 ±0.1 

144 ± 0.4 

0.98T0.00 

— 

45.7 ±0.4 

47.0 ±2.6 

0.97±0.05 


Table 3: The obtained values of anomalous coupling constant a, the mass differ¬ 
ences AM = Mm* — Mm (MeV) of the vector and pseudoscalar mesons, and the 
ratio R = AM{Md)/AM{M,). 




Dd 

m 

d: 

Ds 

Dt 

D., 


Bd 


B* 

Bs 

Bl 

B, 

Scora 

2.85 

2.10 

1.36 

2.67 

1.96 

1.36 

2.32 

2.06 

1.13 

2.15 

1.89 

1.14 

Eichten 

3.07 

2.27 

1.35 

2.90 

2.12 

1.37 

2.60 

2.32 

1.12 

2.42 

2.14 

1.13 

Hagiwara 

3.26 

2.48 

1.31 

3.08 

2.34 

1.32 

2.67 

2.40 

1.11 

2.50 

2.23 

1.12 

Model D 

3.25 

2.47 

1.32 

3.06 

2.35 

1.30 

2.94 

2.65 

1.11 

2.75 

2.47 

1.11 

Average 

3.11 

2.33 

1.34 

2.93 

2.19 

1.34 

2.63 

2.36 

1.12 

2.46 

2.18 

1.13 


Table 4: The sizes of mesons (= 3/(2/5)) in the unit of GeV“^ (1 GeV“^ = 

0.197 fm) obtained from the values of /3’s given in Table 2, and the ratios of the 
vector and pseudoscalar meson sizes. 


10 

































































Figure Captions 


Fig. 1. The meson masses E{j3) = M as functions of the variational parameter (3 
for (a) Dd, (b) Dg, (c) and (d) Bg mesons. 

Fig. 2. The difference of the vector and pseudoscalar meson masses AM = 
Mm* — Mm as functions of the anomalous coupling coefficient a for D^, Dg, B^, 
and Bg mesons. 
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